On the critical temperature of two-channel Kondo model: Reply to |cond-mat /0009283 
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We discuss the comment by Zarand and Zawadowski ( [cond-mat / 0009283 ) on our preprint cond 



mat/0007430 where it has been shown that the strong coupling regime for the two-channel Kondo 



model of two level system can be never realized for any realistic microscopic description. The authors 
of the comment state that the Kondo temperature can be substantially increased due to electron- 
hole asymmetry. Here we show by direct calculation that the electron-hole asymmetry does not 
enter the leading logarithmic approximation. Consequently, we disagree with the aforementioned 
comment. 



Hereby^we discuss the comment by Zarand and Za- 
wadowskiU regarding our preprinto on the Kondo tem- 
perature, Tk, of the two-channel Kondo model based on 
a heavy particle in a double-well potential (DWP) (see 
Ref. HH for a review). A serious deficiency of the original 
2 channel Kondo (2CK) model was too small tunneling 
constants calculated in leading logarithmic approxima- 
tion turns out to be below 0.1 K Authors of Ref. || 
attempted to "cure" the situation by including virtual 
processes via the third excited state in DWP in hope 
that it would increase the tunneling amplitude. 

In Ref. | we demonstrated that when all of the exited 
states are included the Kondo temperature gets reduced 
even further. It turns out that the Kondo logarithmic 
singularity should be cut-off by £3 ~ Tittun <C £f Here 
£3 is the spacing between the second and the third levels 
in DWP, ttun it the tunneling time, and ep denotes the 
Fermi energy. Note that long ago Kagan and Prokof'ev 
made a similar statement in connection with diffusion of 
heavy particles in metals. Since Tk is proportional to the 
cut-off and £3 <C ep, the Kondo temperature turns out to 
be hopelessly small for the strong coupling 2CK limit to 
be achieved experimentally. Moreover, we found that the 
spacing between first two levels always exceeds Tk- This 
prevents realization of this strong coupling regime even 
in principle.Lj In Ref. || we emphasized the importance 
of taking into account all of the exited states as well as 
the fact that the exponential smallness of the tunneling 
probability cannot disappear. 

Formally, this smallness follows from the cancellation 
of virtual contributions of different excited states. At 
the same time, the dramatic reduction of the cut-off and 
thus oITk has a transparent physical reason. Indeed, the 
conventional Kondo Hamiltonian assumes infinitely quick 
spin flip. Therefore, the mapping of the DWP problem to 
the Kondo problem can be valid only for energies below 
fo/ttun - inverse time, which the tunneling process takes. 
Otherwise one could imagine that the Kondo effect takes 
place for a heavy particle in an arbitrary potential, not 
necessarily a double- well one. The reason for the ab- 



sence of the Kondo effect for, e.g., harmonic oscillator 
is that the spacing between the first two levels, A, is of 
the same order as £3 ~ h/t tun . Only the energy interval 
A <C £ <C £3 is available for the formation of the log- 
arithmic singularity! The fact that £F^tun S> Ti follows 
directly from the large mass of the tunneling particle. 
Indeed, the magnitude of DWP should be of the usual 
atomic order ~ 1 V. At the same time, the atomic mass 
exceeds the mass of an electron by about four orders of 
magnitude. Accordingly the tunneling time of an atomic 
particle, t tun should be about two orders of magnitude 
longer than the electron tunneling time scale ~ tijep. 
It is slow tunneling of heavy particles (rather than some 
electron-hole symmetry) that causes the suppression of 
the Kondo temperature. 

The authors of Ref. [j] try to revitalize 2CK model by 
including finite electron-hole asymmetry. Up to our un- 
derstanding, they still do not fully appreciate the fact 
that the proper cut-off energy cannot exceed £3 <C £f 

Below we use an approach, which differs from Ref. |^ 
and show that both electron and hole contributions taken 
separately diverge logarithmically, the ultraviolet cutoff 
being £3 and, therefore, the prefactor in the expression 
for Tk has nothing to do with electron-hole asymmetry. 

The second order corrections to the electron-assisted 
(5A e ) and hole-assisted {5 Ah) tunneling amplitudes can 
be expressed through Matsubara Green functions of elec- 
trons and movable atoms, G e (r) and G a t(r; x\, X2), re- 
spectively: 



5 A e 



l/T 



dr G at (r; -a, a)G e (±r) . 



(f) 



Here we omitted numerical prefactors, put h = 1 and 
used the fact that G e (r) = Gh{— r). We also ignored 
the spatial coordinates and spin indices in G e (r). The 
Green function of the movable atom, G a t(r; X\, x%), is 
determined by the equation 
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d T - U{x) 



G at (r;x,x')=S(T)S(x-x'). (2) 
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We assumed for simplicity that DWP, U (x) , is symmet- 
ric, i.e., U(—x) = U(x), and has the two minima at 
x = ±a. At t < the Green function, G a t(r), van- 
ishes at r < 0, because, according to the model, there 
is only one movable defect. The asymptotic behavior of 
the electron Green function at e ft 3> 1 is (omitting non- 
essential for present discussion constant factors) 



G c i(r) = v{n) — 



ttT dv 
sin ttTt dfi \ s m 



ttT Y 



(3) 



where v(fi) is the density of states at the Fermi level. 
The second term manifests the electron-hole asymmetry, 
whatever it means. Notice, that when r increases, the 
second term decays faster than the first one. Therefore, 
this asymmetry becomes irrelevant at large enough t. 
The r -dependence of the DWP Green function, G at , 
can be found in semiclassical approximation, which is 
valid for excited states as well as for the exponentially 
weak tunneling: 



G at (r) = -B(r)exp {-S c \{xi,x 2 ,t) 



(4) 



where S c \ is the classical action of the particle in the 
inverted DWP: 



S c \{x 1 ,x 2 ,t) 



M fdxV TT . , 



(5) 



5S/5[x(ti)] — 0, x(ti = t) = X2, x(t = 0) = x\. 

One can evaluate the prefactor B(t) in Eq. (||). How- 
ever, we do not need it here. We substitute £1,2 = ±a 
and rewrite Eq.(^) as 



S c i = -tE{t)+ dx y / 2M[E + U(x)} . (6) 

J — a 

Here i?(r)-function is implicitly given by the equation 



dx . 
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(7) 



It follows from Eq. (0), that there are two tunneling 
regimes, with different dependences E(t). The crossover 
takes place at the time scale 



ttnn = I dx 
Jo 
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U(x) - U{a) a - x V U"{a) 



(8) 



The quantity itun ~ hjej, has a meaning of a typical time 
the particle spends under the barrier during the tunnel- 
ing process. Using Eqs. (||) and (Q) one can obtain the 
following asymptotics of S c i at small r 



Sr] = 



2Ma z f a dx __. . 2 , 



(9) 



Equation (9) holds for r -c ^tun- In the opposite case, 



S c i = -tU{o) + / dxy/2M[U(x) - U(a)} 



+2y/ma 4 U"(a) exp 



Jjr/(a) 

M 



(Uun ~ T) 



(10) 



The first term in Eq. (|10|) is nothing but the leading in 
h term in the energy of the ground state which cancels 
out for the on-shcll processes. The second term is the ac- 
tion corresponding to the tunneling from one minimum to 
the other. The classical trajectory for this process corre- 
sponds to the kink with the characteristic time scale itun- 
The last term in the action characterizes the corrections 
due to the finite size of the kink (we write it here only 
for the references purpose). Substituting Eqs. (||) and 
Eq. (|) into Eq. (g), and using the condition e^tim 3> 1, 
and asymptotic expansions (||) and ( |l0| ) we obtain 
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dx^2M[U(x) - U{a)} 
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0(1) . (11) 



We are not discussing here if the two contributions, SA e 
and 6 Ah, add or cancel each other (partially or com- 
pletely). What we would like to emphasize is that for 
each of the two contributions the cut-off of the logarith- 
mic dependence of each contribution is determined only 
by the tunneling time. This time exceeds inverse Fermi 
energy of the electrons by orders of magnitude, since the 
atomic mass M is much bigger than the electron mass. 

In fact, this asymmetry (i.e. finite dv jd\i) is relevant 
only at short times r -C itun, where the integral over r in 
Eq. |l]) is not logarithmic. We conclude that the origin 
of ultraviolet limit of the logarithm has nothing to do 
with the electron hole asymmetry. The smallness of the 
prefactor in the expression for the Kondo temperature 
expresses the simple fact that for the logarithmic diver- 
gence to appear in any physical property a heavy particle 
must tunnel and its tunneling time is finite, just as it was 
explained in Ref. ||. At times smaller than the character- 
istic time of tunneling, the heavy particle does not feel 
the two- well structure of the potential yet, and does not 
acquire any logarithmically divergent correction. Here 
we demonstrated this fact by an explicit calculation al- 
ternative to Ref. [| The physical model used in Ref. [j] 
does not differ from and is the same as the one that was 
discussed above, of Ref. |l|. For these reasons we believe 
that the calculation is erroneous. 

The final remark of Ref. |l| that the degeneracy of the 
low-lying states may be "guaranteed by some local sym- 
metry (rotational or possibly time reversal)" is also, in 
our opinion, irrelevant for metallic systems. Time re- 
versal symmetry is known to guard the double degener- 
acy only for systems, which Hamiltonian contains a half- 
integer spin (Kramers degeneracy). Such a Hamiltonian 
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leads to the regular Kondo problem and has nothing to 
do with the original two-level system. As to rotational 
symmetries, they cause no degeneracies for the ground 
state, because the ground state of a single atom always 
belongs to the one- dimensional representation A\ or A 
of the point symmetry group. 

To summarize, we disagree with the statement of 
Ref. 0, that the suppression of the Kondo temperature 
is a consequence of the electron hole symmetry. We in- 
sist that the Kondo temperature for a particle in DWP is 
always parametrically less than the splitting between the 
two lowest state. Accordingly, the strong coupling 2CK 
regime never takes place in metals. 
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